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Spatially extended coherence induced by pairing correlation 
in low-frequency vibrational excitations of neutron drip line nuclei 
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Heavy Ion Nuclear Physics Laboratory, RIKEN, Hirosawa 2-1, Wako, Saitama 351-0198, Japan 

(Dated: February 9, 2008) 

Role of pairing correlation for emergence of low-frequency vibrational excitations in neutron drip 
line nuclei is discussed paying special attention to neutrons with small orbital angular momentum I. 
Self-consistent pairing correlation in the Hartree-Fock-Bogoliubov (HFB) theory causes the change 
of the spatial structure of the quasiparticle wave functions; "the pairing anti-halo effect" in the 
lower component and "the broadening effect" in the upper component. The resultant spatial distri- 
bution of the two-quasiparticle states among \ow-l neutrons, "the broad localization", brings about 
qualitatively new aspects, especially the large transition strength of the low-frequency vibrational 
excitations in nuclei close to the neutron drip line. By performing HFB plus quasiparticle random 
phase approximation (QRPA) calculation for the first 2 + states in neutron rich Ni isotopes, the 
unique role of self-consistent pairing correlation is pointed out. 



PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Study of low-frequency vibrational excitations in neu- 
tron drip line nuclei is one of the most interesting sub- 
jects in nuclear structure physics. Vibrational excita- 
tions are microscopically represented by coherent super- 
position of two-quasiparticle states (or one-particle - one- 
hole states in closed shell nuclei) In stable 
nuclei, the main configurations of the low-frequency vi- 
brational excitations are two-quasiparticle states among 
tightly-bound single-particle states. In neutron drip line 
region, by contrast, the contributing single-particle states 
are tightly-bound states, loosely-bound states, resonant 
and non-resonant continuum states. Consequently the 
two-quasiparticle states among them have rich variety of 
the spatial structure. Therefore the low-frequency vibra- 
tional excitations, as the coherent motions among such 
two-quasiparticle states, may have qualitatively new as- 
pects compared to those in stable nuclei. 

We expect unique impacts of the two-quasiparticle 
states involving loosely-bound low-£ neutrons on the low- 
lying excitations. Loosely-bound \ow-t neutrons have 
an appreciable probability to be outside the core nu- 
cleus that leads to the neutron halo structure (see Refs 
Hi [1 LH LI- El f° r reviews), and the coupling to the nearby 
continuum states causes the soft excitation with large 
transition strength as a non-resonant single-particle ex- 
citation • An example is the soft dipole excitation in 
light halo nuclei 0, S (13, EH E3, E3, E3, EE EE 113 , typ- 
ically in n Be and 11 Li, where significant El strength is 
observed. 

On the other hand, the role of two-quasiparticle states 
among loosely-bound low-£ neutrons and resonant states 
for low-frequency vibrational excitations is not well un- 
derstood so far. It is still an open question whether the 



coherent motions among such two-quasiparticle states 
can be realized in spite of the variety of the spatial struc- 
ture, and the resultant low-lying excitations may possess 
novel features in neutron drip nuclei. 

In the present study, we emphasize the unique role of 
self-consistent pairing correlation that changes the spa- 
tial structure of the quasiparticle wave functions; "the 
pairing anti-halo effect" in the lower component and "the 
broadening effect" in the upper component. We show 
that the resultant two-quasiparticle states can cause the 
coherence among them not only in the surface region but 
also in the spatially extended region. Consequently the 
spatially extended coherence induced by pairing correla- 
tion leads to the large collectivity of the low-frequency 
vibrational excitations. 

In Sec|nJ by solving the Woods-Saxon potential plus 
HFB pairing model, we examine the pairing correla- 
tion in loosely bound nuclei and the induced change 
of the spatial structure of the quasiparticle wave func- 
tions. In Sec lIIII we examine the spatial structure of 
two-quasiparticle states in neutron drip line nuclei. In 
Sec lIVI we perform Skyrme-HFB plus QRPA calculation 
for the first 2 + states in neutron rich Ni isotopes. By 
comparing three types of calculations; HFB plus QRPA, 
resonant BCS plus QRPA, and RPA, we emphasize the 
unique role of self-consistent pairing correlation for real- 
izing low-frequency vibrational excitations and the large 
collectivity in neutron drip line nuclei. Conclusions are 
drawn in SeclVl 



II. SPATIAL STRUCTURE OF 
QUASIPARTICLE WAVE FUNCTIONS 



Model 
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For our qualitative discussion, we solve the HFB equa- 
tion in coordinate space 

EI El HI EH with spherical 
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symmetry, 

-,(h l3 {ry) A(r,r') \( Ulj {E,r>) 
\ A(r,r') -htjiry) J^v^Ey) 

E + \ \ f ui 3 (E,r) 
E-X J yv h (E,r) 

where E is the quasiparticle energy, A is the Fermi energy, 
and uij(E,r) (vij(E,r)) is the upper (lower) component 
of the radial quasiparticle wave function. The upper and 
lower components are the generalization of the uij and 
vij coefficients in the BCS approximation. 

For the mean-field hamiltonian hij(r,r'), the Woods- 
Saxon pote ntial together with the related spin-orbit po- 
tential (H, 



A 



(a-l), (2) 



where 



f{r) = 1/ 1 + exp 



Rws 



(3) 



is adopted. A = h/mc is the reduced Compton wave- 
length of nucleoli. The parameters are fixed to be a = 
0.67 fm and v — 32 in accordance with Ref.[22j. The 
strength Vws and the radius Rws are chosen to simulate 
the shell structure of nuclei under consideration. For the 
pairing correlation, we self-consistently derive the pairing 
potential A(r, r') from the density dependent zero-range 
pairing interaction, 
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that leads to the local pairing potential 
A(r,r') = A(r)S(r-r'), 

with 
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where P a is the spin exchange operator. The normal 
density p(r) and the abnormal density p(r) are defined 
by 

p(r) = iE( 2 J + 1 )K(^».'-)} 2 . ( 7 ) 
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To measure the pairing correlation in nuclei under con- 
sideration, we calculate the pairing energy, 



E pair = - J A{r)p{r)^r 2 dr. 



(9) 



For comparison with Ref. 22J,|23J,|24|, the average pairing 
gap defined by 

/ A{r)p(r)r 2 dr 



A = 

is also calculated. 



J p{r)r 2 dr 



(10) 



> 



> 

s 



— • 


• — 


Self-consistent pairing potential 




F =-680 MeVfm" 3 
















y 



FIG. 1: The average pairing gap A and the pairing energy 
E pa i r obtained by the Woods-Saxon potential plus HFB pair- 
ing model are plotted as a function of \=E3 3l/2 . The density- 
dependent pairing force with V pa ir~— 680 MeV fm -3 is used. 
The box size is fixed to be 75 fm. 



B. Inputs 

We examine the spatial structure of the two- 
quasiparticle states of neutrons around 86 Ni as the illus- 
trative example of neutron drip line nuclei. The nucleus 
86 Ni is at the neutron drip line within Hartree-Fock (HF) 
calculation with Skyrme SLy4 force 25], and has the 
loosely-bound 3si/2 state as the Fermi level (see Fig ll6H . 
The Woods-Saxon potential with Vws — —41.4 MeV and 
Rws = 5-5 fm simulates the neutron shell structure, and 
the obtained single-particle states around the Fermi level 
are lg 9/2 (e h = -4.93 MeV), 2d 5/2 (e h = -1.77 MeV), 
3si/ 2 (e h = -0.50 MeV), resonant d 3/2 (e p « 0.14 MeV), 
and resonant g 7 / 2 (s p « 0.98 MeV). The index p (h) rep- 
resents all necessary quantum numbers to specify the par- 
ticle (hole) state. In the present study, the box boundary 
condition with a box size Rbox — 75 fm is imposed. We 
elaborately examine the applicability of the box bound- 
ary condition for describing pairing correlations in loosely 
bound nuclei in the following subsections. 

We solve Eq.QJ only for neutrons, and p(r) (p(r)) rep- 
resents the neutron normal (abnormal) density in the 
present section. The pairing strength is fixed to be 
V pair = -680 MeV fm" 3 with cut-off energy E cut = 50 
MeV. The parameter p c is taken to be 0.08 fm -3 (« p(0)) 
for the surface-type pairing field. The obtained average 
pairing gap is A « 1.3(» 12/^86) MeV. 



C. Pairing correlation in loosely-bound nuclei 

In the present subsection, we discuss the pairing cor- 
relation in loosely-bound nuclei. We focus on the ap- 
plicability of HFB calculation with the box boundary 
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FIG. 2: The box size dependence of the average rms ra- 
dius of the lower component of the S1/2 quasiparticle wave 
function and the pairing energy obtained by the Woods- 
Saxon potential plus HFB pairing model with the condition 
of A = e 3si = -0.052 MeV. 
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FIG. 4: The occupation probability of the S\n state with the 
condition of A=£3s 1/2 =— 0.052 MeV obtained by the Woods- 
Saxon potential plus HFB pairing model. The results ob- 
tained with Rtox = 75 and 400 fm are compared. 
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FIG. 3: The average effective pairing gap of the S1/2 state as a 
function of A = £3s 1/2 • The results of the Woods-Saxon poten- 
tial plus HFB pairing model with the self-consistent pairing 
potential and the analytic pairing potentials of A = 1.0 and 
1.5 MeV are compared. 



condition for the extreme situation; the pairing effect 
on the S1/2 state in the limit of e Sl/2 — > 0. In FigQ] 
the average pairing gap and the pairing energy around 
86 Ni are plotted as a function of A = £3s 1/2 ■ The single- 
particle energy changes by varying Vws while keeping 
the other parameters fixed. The box boundary condi- 
tion with Rbox = 75 fm is imposed. Pairing correlation 
becomes monotonically stronger as £3s 1/2 increasing up 
to around -0.5 MeV, because of the change of the shell 
structure that is self-consistently taken into account in 
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FIG. 5: The abnormal density p[r) and the contributions 
from the low-lying resonant states pij(r) obtained by the 
Woods-Saxon potential plus HFB pairing model in the case 
of A=£3s 1/2 =— 0.052 MeV. The box size is fixed to be 75 fm. 



the pairing potential. The pairing correlation is much 
enhanced in the very loosely-bound region, A = £3 Sl/2 > 
—0.5 MeV, due to the increasing coupling to the nearby 
continuum 26]. The importance of the coupling to the 
non-resonant continuum states for pairing correlation is 
also discussed in connection with di-neutron correlation 

To examine whether such enhancement of pairing cor- 
relation is due to the artifact of the box boundary con- 
dition, we examined the box size dependence up to 
Rbox — 400 fm. In Fig|2]the pairing energy in the case of 
A = £3s 1/2 = —0.052 MeV is shown as a function of Rbox- 
The strength Vws is fixed to be -39.0 MeV. We confirm 
that the box size 75 fm is enough to describe the pairing 
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FIG. 6: The box size dependence of ps 1/2 (r) in the case of 



-0.052 MeV. The results obtained with R box 



30, 75, and 400 fm are compared. 



correlation in such very loosely-bound nuclei. 

To estimate the pairing effect in the 3si/2 state, 
Hamamoto et al. [23 introduced the effective pairing gap 
by A e //(s 1 / 2 ) = E Sl/2 , where the smallest quasiparticle 
energy E Sl/2 is calculated with the condition A = £3s 1/2 
for the discrete state satisfying (A + E Sl/2 ) < 0. In the 
case of (A + E Sl/2 ) > 0, the value of E Sl/2 is obtained by 
calculating the derivative of the phase sift of the quasi- 
particle wave function u Sl/2 (E Sl/2 , r). In our analysis, to 
estimate the average pairing effect on the resonant state, 
we introduce the average effective pairing gap by 



Ae//(siy 2 ) = E av (s 1/2 ), 



(11) 



where E av (si/ 2 ) is the average energy of the resonant 
state with the occupation probability, 



Eav{si/2) 



1 



N res (si/ 2 ) 

x J2V(E WI < W)(v Sl/2 , n ) 2 E Sl/2 , n , (12) 

n 

where 9(x < a) =1 for x < a, otherwise 0. The occupa- 
tion probability at Eij_ n is defined by 

(%-.„) 2 = / {vy(E Ii>n ,r)} 2 dr. (13) 
Jo 

The integrated occupation probability over the resonant 
state, 



N^=Y l 0(E lj ,n<W)( 



Vlj,n) 2 , 



(14) 



is introduced. 

In Fig[3| A e ff(si/2) is plotted as a function of A = 
£3s 1/3 - As A = £3s 1/2 increases, A e ff{si/ 2 ) monotonically 
increases. The enhancement of A e //(si/2) in the limit of 



A = £3s 1/2 — > appears. This behavior is consistent with 
the pairing correlation of the total system measured by 
A and E pa ir (see FigQ]). This means that the strong 
coupling between the pairing potential and the Si/ 2 state 
remains in the limit of A = £3s 1/2 — > 0. 

We also examined the box size dependence of 
A e //(s 1 / 2 ). In Fig 0| the occupation probability of the 
Si/2 state in the case of A = £3 Sl/2 = —0.052 MeV 
is shown as a function of E Sl/2 . The results obtained 
with Rbox — 75 and 400 fm are compared. Although the 
level density of the discretized continuum states is very 
different, the average structure of the resonance is un- 
changed. Namely the integrated occupation probability 
N, 



(res) 

1/2 



is 0.430 (0.429), and the average effective pairing 
gap A e// (s 1/2 ) is 1.15 (1.13) MeV with R box = 75 (400) 
fm. And also, as we will discuss in the next subsection, 
the average root-mean-square (rms) radius of the lower 
component of the s 1 / 2 state is independent of Rbox larger 
than 50 fm (see FigEJ. 

To emphasize the importance of self-consistency in the 
pairing potential, we perform HFB calculation with the 
analytic surface- type pairing potential |22j . 



A(r) oc r 



df(r) 
dr 



(15) 



where f(r) is the Woods-Saxon type function of Eq.@. 
The average pairing gap with the analytic potential is 
defined by 



A = 



g*™ A(r)f(r)r 2 dr 



(16) 



In Fig0the average effective pairing gap A e //(s 1 / 2 ) 
obtained by HFB with the analytic pairing potentials of 
A = 1.0 and 1.5 MeV are also shown. These A e //(s 1 / 2 ) 
decrease rapidly in the limit of zero binding energy, due 
to the decoupling between the Si/ 2 state and the pairing 
potential as discussed in Ref . |24| . contrary to the HFB 
calculation with the self-consistent pairing potential. 

The low- lying si/ 2 state is the dominant component 
of p(r) outside the surface region in the HFB calculation 
with the self-consistent pairing potential. To see this, 
we define the contributions from the low-lying resonant 
states to p(r) by 



Pij(r) 



(2j + 1) 
47rr 2 

n 

xuij(Ei j<n ,r)vij(Eij <n , r). 



0{Eij,n < 10) 



(17) 



In FigO Pij(f) are shown in the case of A = £3si /2 = 
—0.052 MeV with Rbox — 75 fm. As clearly seen, the 
Psi/2 ( r ) h &s a sizable component outside the surface re- 
gion, r > 10 fm. In Fig[f)| the box size dependence of 
Ps 1/2 {r) is examined, and we confirmed that R box = 75 
fm is enough to describe the pairing correlations in such 
very loosely-bound nuclei. Because of the spatially ex- 
tended component in the self-consistent pairing poten- 
tial, the pairing field can act on the Sx/a state, although 
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s lower component 

Self-consistent pairing 

Analytic pairing (1.5 MeV) 

Analytic pairing (1 .0 MeV) 

No pairing 



-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



TABLE I: The quasiparticle energies of the low-lying reso- 
nant states E i f. es ' > obtained by the Woods-Saxon potential 
plus HFB pairing model with Vws = —41.4 MeV. The aver- 
age rms radii of the lower components Rij and the integrated 
occupation probabilities N^ ea ' > and iVy(l.O) are shown. The 
single-particle energies of the corresponding bound or reso- 
nant states Eij are also listed. The box size is fixed to be 75 
fm. 
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FIG. 7: The average rms radius of the lower component of 
the S\/2 quasiparticle wave function as a function of A=£3s 1/2 
obtained by the Woods-Saxon potential plus HFB pairing 
model. The results with the self-consistent pairing potential 
and the analytic pairing potentials of A = 1.0 and 1.5 MeV 
are compared. The radius without pairing is also shown. The 
box size is fixed to be 75 fm. 

the loosely-bound Si/ 2 state is spatially extended. The 
coupling between the self-consistent pairing potential and 
the loosely-bound s±/2 state becomes much stronger by 
repeating the process self-consistently because of the 
pairing anti-halo effect that we will discuss in the next 
subsection. On the other hand, this process is absent in 
the analytic pairing potential, and the coupling to the 
loosely-bound Siy 2 state becomes artificially weaker for 

£ 3s 1/2 0. 

The enhanced pairing energy in the very loosely-bound 
nuclei originates from the spatially extended component 
of p(r). As shown in Fig[21the pairing energy in the case 
of A = e 3si/2 = -0.052 MeV obtained with R box = 400 
(30) fm is -20.368 (-20.238) MeV. The difference 130 keV 
is due to the spatially extended component of p(r) shown 
in FigEl 

D. Pairing anti-halo effect in lower components 

Self-consistent pairing correlation in HFB has a unique 
role that changes the spatial structure of the quasiparticle 
wave functions in loosely-bound nuclei. In the present 
subsection, we focus on the spatial structure of the lower 
component, and the structure of the upper component 
will be discussed in the next subsection. 

According to the asymptotic HFB equation with 
Vij (r, r') — > and A(r, r') — » in the limit of r — > oo, the 
lower component of the quasiparticle wave function de- 
cays exponentially for any quasiparticle energy £y > 0, 

vij{Eij,r) -► exp(-/3yr), (18) 

where /3y = ^2m(Eij — X)/h 2 . The HFB quasiparticle 
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FIG. 8: The effective hole wave functions of the S1/2, ^3/2, 
d 5 / 2 and g 7 / 2 resonant states obtained by the Woods-Saxon 
potential plus HFB pairing model are shown. The functions 
in the resonant BCS approximation are also plotted. The 
strength Vws = —41.4 MeV and the box size Rbox = 75 fm 
are fixed. 

energy By is well approximated by the canonical quasi- 
particle energy, E tj ~ Eff n = ^(e" 1 " - A) 2 + (A™™) 2 , 

where ef"™ and A™" are the canonical single-particle en- 
ergy and the canonical pairing gap of the corresponding 
state plj . This means that /3y stays at finite value 

Aj -> ^2mE h /h 2 > ^2mA™/S 2 > (19) 
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with non-zero Ay on for A — > 0. Therefore self-consistent 
pairing correlation changes the asymptotic behavior of 
vij (Eij , r) , and especially acts against the development of 
the infinite rms radius that characterizes the wave func- 
tions of s and p states in the limit of vanishing binding 
energy. This change of the spatial structure in the lower 
component of the quasiparticle wave function is called 
"the pairing anti-halo effect" in Ref. 28]. 

In Fig0the average rms radii of the lower component 
of the S]/2 state with the self-consistent pairing potential 
and the analytic pairing potentials of A=1.0 and 1.5 MeV 
are shown as a function of A = e^ Sl , 2 . The average rms 
radius Rij is defined by 

y « 

x / 60i r 2 K(Sy.„,r)} 2 dr. (20) 
Jo 

The radius of the 3si/2 state without pairing goes to in- 
finity for £3s 1/2 — > 0. On the other hand, the divergence 
is suppressed in HFB. For A = e^ Sl/2 < —1.5 MeV, irre- 
spective of the pairing potential, the radii obtained by the 
HFB calculations have a similar behavior. As approach- 
ing the zero binding energy, the spatial structure of the 
lower component becomes sensitive to the treatment of 
pairing correlation in accordance with the behavior of the 
average effective pairing gap A e ^/(s!/ 2 ) shown in FigOU 
In connection with the spatial structure of two- 
quasiparticle states, we examine the effective hole wave 
function of the transition with multipolarity L defined by 

/SW=r%-(^>,r-)/%>. (21) 

The factor uy, n (> 0) is introduced for the normalization. 
The effective hole wave functions of the S1/2, £?3/2i ^5/2 
and #7/2 resonant states with L = 2 are plotted in Fig(8] 
The strength Vws is fixed to be -41.4 MeV to simulate 
the neutron shell structure in 86 Ni. The properties of 
these resonant states are summarized in Tabled In HFB 
the effective hole wave function is localized according to 
Eq. (|18|l . irrespective of the corresponding single-particle 

energy. However the localization depends on I. The func- 

(2) 

tion fs 1/2 (r) has the large component not only around 
the surface region, but also the extended region, 10 fm 
< r < 20 fm. The component in the extended region re- 
mains, but becomes smaller in the d 3 /2 an d ^5/2 states, 
and negligible in the g?/ 2 state. In Subsec lTTll we will 
discuss the impact of such spatially extended part of the 
lower component in connection with the new aspects of 
low-frequency vibrational excitations in neutron drip line 
nuclei. 

For comparison, we analyze the effective hole wave 
functions in the resonant BCS approximation. Although 
the BCS approximation is not well defined for unstable 
nuclei due to the unphysical neutron gas [l^, recently 
the extended version, the resonant BCS approximation, 
was proposed [H |H EJ HI • In this approach, resonant 



states are taken into account for the particle states. Ap- 
plicability of the resonant BCS method for ground state 
properties is discussed in Ref. [33J. In this approxima- 
tion, the upper and lower components of the quasiparticle 
wave function are simply proportional to the correspond- 
ing single-particle wave function ¥>y(£y, nj r ), 

uij (Eij,n,r) = u^ c n s tpij (sij, n ,r), , 22 . 
% (Eij,„,r) = v% c n s ipij (£«,„, r). 

The quasiparticle energy £?y, n is related to the single- 
particle energy ey >n by 

Eij,n = yW>- A) 2 + (A z „ „) 2 , (23) 

where the Fermi energy A and the pairing gap Ay jn . 
v Vj C n S i u f^rf) 1S tne occupation (unoccupation) ampli- 
tude. The index n runs over not only the bound states 
but also the resonant states. For r — + 00, the radial 
single-particle wave function of the bound state decays 
exponentially, 

<Plj(eij,n,r) exp(-ay,„r), (24) 

where ay. rl = ■J— 2mei^ n jW- . On the other hand, the 
wave function in the continuum region satisfies the scat- 
tering boundary condition, 

ipy(eij,r) -> [cos(5ij)rji(kijr) - sin(5y)rnj(fcy70i5a) 
- sin(%r + <5 y ) (25b) 

where fcy = y/2msij/h 2 . ji and ni are spherical Bessel 
and Neumann functions, and Sij {5[- = 6ij — In/ 2) is the 
phase shift corresponding to the angular momentum (Ij). 

As shown in FigJSl the effective hole wave functions of 
the d 3 / 2 and 37/2 resonant states in the resonant BCS 
approximation diverge for r — ► 00 in contrast with those 
in HFB. The effective hole wave functions of the bound 
single-particle states in the (resonant) BCS approxima- 
tion are essentially the same with those without pairing. 

(2) 

The function for the halo state, fs 1/2 (r), is spatially very 
extended due to the absence of the pairing anti-halo ef- 
feet. On the other hand, fd s/2 ( r ) nas almost the same 
spatial structure with that in HFB, because the corre- 
sponding single-particle state is the bound state without 
halo. 



E. Broadening effect in upper components 

Self-consistent pairing correlation also changes the spa- 
tial structure of the upper component of the quasiparticle 
wave function. Asymptotic behavior of the upper com- 
ponent depends on the quasiparticle energy |18|. In the 
discrete region (£^y, n + A) < 0, the wave function decays 
exponentially, 

uij(Eij in ,r) -> exp(-ay )n r), (26) 
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FIG. 9: The localization indicator of the resonant d 3 / 2 state 
at £ (res) = 1.25 MeV in the Woods-Saxon potential plus HFB 
pairing model is shown. The results obtained with Rbox = 30, 
75, 150, and 400 fm are compared. 




g State HFB E = 2.05 + 0.5 MeV 

"BCS" 8= 0.98 + 0.5 MeV 

State HFB 1-25 ±0.5 MeV 

- "BCS" s= 0.14 ± 0.5 MeV 




FIG. 10: The localization indicators of the resonant state 
at E {res) = 1.25 MeV and the resonant g 7/2 state at E {rea) = 
2.05 MeV in the Woods-Saxon potential plus HFB pairing 
model are shown. The "BCS" localization indicators of the 
corresponding states are also plotted. The box size is fixed to 
be 75 fm. 



where OLij, n = — 2m{Eij,n + X)/h 2 . In the continuum 
region (Ey + A) > 0, on the other hand, 



uij(Eij,r) 



[cos (Sij)rji(kijr) 
sin(%r + 6'fj), 



sin (Sij) rni(kijifflp) 
(27b) 



with ky = y/2m(Eij + \)/h 2 . 

In the neutron drip line region with A — > 0, all quasi- 
particle states are in the continuum region. To examine 
the spatial structure of the upper components, we intro- 
duce the localization indicator, 



0.15 
0.12 
0.09 
0.06 
0.03 
0.00 



* state 



HFB E = 0.80 ± 0.5 MeV 
"BCS" i'=-0.50McV 




HFB E= 1.72 ±0.5 MeV 
"BCS" e= -1.77 MeV 




FIG. 11: The same as Fig llOl but for the resonant S1/2 
MeV. 



0.80 MeV and the resonant d 5/2 state at E {res) 



state at 
1.72 



The summation is taken around the resonance energy 
E^ r - es ^ within the energy window SE. In principle, the 
localization indicator is independent of Rb ox with arbi- 
trary 5E, if Rbox is enough large. In FigOHthe indicator 
of the d 3/2 state at E^ = 1.25 MeV obtained with 
Rbox — 30, 75, 150, and 400 fm are compared. The en- 
ergy window 5E is fixed to be 1.0 MeV. Apart from the 
artifact around the box boundary, these results coincide 
in good accuracy. The box size Rbox — 75 fm is enough 
for our discussion on the spatial structure of the two- 
quasiparticle states around 86 Ni, because it is required 
to describe the spatial structure of the upper components 
within 30 fm as seen in Fig|S] 

In Figmi the indicators of the d 3 / 2 state at 
E (res) =1 25 MeV and the g 7/2 state at E (res ^ =2.05 
MeV, that originate from the resonant single-particle 
states, are shown. The occupation probabilities within 
the energy window, 



\E, 



lj,n 



E, 



{res) I 

y 



<SE/2)(v^ n ) 2 , (29) 



together with SE = 1 MeV are listed in Tabled To ex- 
amine the role of self-consistent pairing correlation, as 
the reference distribution, the localization indicator in 
the resonant BCS approximation (we call the "BCS" lo- 
calization indicator) for e,' es ^ > is introduced by 



Lij(r;SE) = 



(l-JV w (tf£)) x 

-4 es) |<<5i?/2){^(^,„,r)} 2 .(30) 



L l3 (r;6E) 



£< 



\E, 



lj,n 



E/r'\<6E/2){u lj (E lj , n ,r)}\28) a 



The single-particle wave function tpij(sij,ni r ) with ey jn > 
is the solution of the Woods-Saxon potential of Eq.J2J) 
with the box boundary condition. 

In the <77/ 2 resonant state, the difference between the 
FB and "BCS" indicators is small within 20 fm. On the 
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other hand, in the state, the spatial localization of 
the resonance is weakened in HFB; namely, the compo- 
nent within the centrifugal barrier (r < 10 fm) is reduced, 
and enhanced in the extended region (10 fm < r < 20 
fm) . This is the direct manifestation of the strong mixing 
with the nearby non-resonant continuum states. As dis- 
cussed in Ref. [l8l I22I l33l l34j , the resonance width of the 
quasiparticle state originating from the resonant (bound) 
single-particle state is broadened (acquired) in HFB. This 
effect is more prominent in \owev-l states with smaller 
quasiparticle energy. Consequently, the spatial structure 
of the upper component of the quasiparticle wave func- 
tion changes, and the spatial localization is weakened by 
the mixing with continuum states. To emphasize this 
effect induced by self-consistent pairing correlation; the 
broadened resonance width and the spatially weakened 
localization, we call "the broadening effect" in the upper 
component. 

In FigEJthc localization indicators of the S1/2 state 
at E (res 1 = 0.80 MeV and the d 5/2 state at E ires 1 = 
1.72 MeV, that originate from the bound single-particle 
states, are shown. The "BCS" localization indicator with 
E lj,n < is defined by 

L«(r) = (1 - N lj {5E)){ Vlj {ei j , n ,r)f. (31) 

In these states, the spatial structure of the HFB and 
"BCS" indicators is very different. The "BCS" indicator 
is localized, but non- localized in HFB, because of the dif- 
ferent boundary conditions, and the broadening effect is 
absent in the (resonant) BCS approximation. The broad- 
ening effect is strong in the upper component of the S1/2 
and c?3/2 states. Especially the localization around the 
surface region is strongly affected in the S1/2 state. 

III. SPATIAL STRUCTURE OF 
TWO-QUASIPARTICLE STATES 

A. Two-quasiparticle states in HFB 

The spatial localization of two-quasiparticle states (lp- 
lh states in closed shell nuclei) is a necessary condition 
to realize the coherent motions among them. In stable 
nuclei the spatial localization of the contributing two- 
quasiparticle states is realized around the surface region 
(for example, see Ref.|35l l36| ). In neutron drip line nu- 
clei, on the other hand, because the contributing single- 
particle states are tightly-bound states, loosely-bound 
states, resonant and non-resonant continuum states, the 
spatial distributions of the two-quasiparticle states have 
rich variety of the spatial character. Therefore, first of 
all, the realization of coherent motions among them is not 
obvious. If they can be realized, we may expect quali- 
tatively new aspects of low-frequency vibrational excita- 
tions in neutron drip line nuclei. 

In the present study, we analyze the transitions from 
the ground state to excited states within the same nu- 



cleus. The spatial distribution of the particle-hole com- 
ponent of the two-quasiparticle state between the up- 
per component uy y (Ey y > r) and the lower component 
vij (Eij , r) is defined by 

FivHr) = u Vi >(E Vj ,,r)r I 'vi j {E lj ,T)/N wl , (32) 

where N uv — yj\ — (vi'j<) 2 vij (> 0) is the normalization 
factor. 

One of the novel features of low-frequency vibra- 
tional excitations in neutron drip line nuclei is that two- 
quasiparticle states composed of quasiparticles both in 
the continuum region are inevitably involved. In HFB, 
owing to the pairing anti-halo effect, the asymptotic be- 
havior for r — ► 00 is, 

F&Hr) - r L sm(k VJ ,r + S' Vi ,) exp(-/?yr). (33) 

Since /3y stay at finite with non-zero pairing correla- 
tion, the spatial distribution is localized irrespective of 
the quasiparticle energy. 

Recently QRPA calculation with the resonant BCS ap- 
proximation is applied to describe vibrational excitations 
in loosely-bound nuclei [8Sl E^j. In this method, ac- 
cording to Ea. (|25b|) . the spatial distribution of the two- 
quasiparticle state composed of quasiparticles both in the 
continuum is non-localized function with the asymptotic 
behavior, 

F&\r) - r L ^{k vf r + 5 l vr )s\n(k l3 r + 5[ 3 ). (34) 

Consequently the spatial integration, namely the transi- 
tion matrix element, diverges. In practice this divergence 
may be suppressed by imposing the box boundary con- 
dition with small box size, typically Rbox = 15 fm |38| . 
However a large box radius is needed to represent the 
wave functions with the spatially extended structure in 
good accuracy. Therefore this approximation is uncon- 
trollable for the description of excitations in nuclei close 
to the neutron drip line. In Ref. 39] the low-lying dipolc 
modes in 26,28 Ne are studied by means of quasiparticle 
relativistic RPA with the resonant BCS approximation. 
In this calculation, although the resonant states are cal- 
culated by imposing the scattering boundary condition, 
the obtained results don't suffer with the problem of the 
divergence because of the peculiarity of the shell struc- 
ture. For example, the two-quasiparticle state between 
the resonant ff/ 2 an d P3/2 states doesn't contribute to 
the dipole excitation in 28 Ne. Moreover the unoccupied 
states outside of the pairing active space are obtained by 
expanding on a set of the harmonic oscillator basis, and 
the divergence of the two-quasiparticle states between the 
occupied part of the resonant states and the non-resonant 
continuum states is artificially eliminated. 

B. Broad localization of two-quasiparticle states 

As the illustrative example of neutron drip line nu- 
clei, we examine the spatial structure of the neutron 
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FIG. 12: The (r) between the bound states and the res- 
onant states with e p — £h<5 MeV in 86 Ni. 



FIG. 14: The F„„ (r) associated with the low-lying Si/ 2 , ^3/21 
d 5/ / 2 , and g 7 / 2 resonant states in 86 Ni obtained by the Woods- 
Saxon potential plus HFB pairing model. 




FIG. 13: The first moments R$ of F$ (r) in FigGU The 



regions \R$ 



(2) 



quadrupole two-quasiparticle states in 86 Ni. To charac- 
terize the spatial structure of Fuv (r), we introduce the 
first moment, 



(r)|o?r 



f 
Jo 



(r)\dr 



and the variance, 



f Rb °*r*\F^(r)\dr 



(L) 



r 

Jo 



(r)\dr 



\ 1 L uv J 



(35) 



(36) 



If the regions \R 



(L) 



r\ < a^uv 



have overlap, we expect 




/ / / / / / / / / / 

^ ^ ^ ^ & ^ ^ ^ # i 

FIG. 15: The first moments 7?S of F,$(r) in FiglTl The 
regions \Ruv — r\ < a^v are shown by vertical line segments. 
The result without the broadening effect in the upper com- 
ponents of the quasiparticle wave functions is also shown. 



(a) 3s 1/2 



resonant d 



3/2; 



(b) 2d 



5/2 



resonant d 



and (c) 2^5/2 — > resonant 37/2. The corresponding R 



3/2i 
(2) 



and the regions \R. 



(2) 

ph 



< a 



(2) 
ph 



ph 



are shown in Fig El 



the correlations among the two-quasiparticle states. 

In Fig^| the spatial distributions without pairing, 
Fph ( r ) ' are shown- Because the available particle states 
are only the d 3 / 2 an d 37/2 resonant states, the lp-lh ex- 
citations to the resonant states with e p — Eh < 5 MeV are localized in HFB. The corresponding Ruv and the regions 



These configurations have overlap around the surface re- 
gion, and the correlations among them bring some low- 
lying collectivity. However the lowest lp-lh state (a) only 
have the appreciable component outside the surface re- 
gion that can't correlate with the other configurations, 
and the lowest RPA solution becomes a single-particle 
like excitation as we will discuss in Subsec llVBI 

In Figll4l the spatial distribution of the two- 
quasiparticle states associated with the low-lying S1/2, 
^3/2 1 ^5/2 an d 37/2 resonant states are plotted. Due to 
the pairing anti-halo effect, all distribution functions are 

(2) 
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|p(2) 



r\ < u^J are shown in Fig ll5l The distributions 
of the two-quasiparticle states between the s 1 / 2 , d 3 / 2 an d 
d-5/2 resonant states have the sizable components in the 
spatially extended region around 10 fm < r < 15 fin, 
in addition to the surface region where the localization 
is achieved in stable nuclei. Because of the spatially ex- 
tended structure, the transition matrix elements can be 
larger than those of the tightly bound states in stable 
nuclei. Therefore we can expect the large collectivity 
of low-frequency vibrational excitations by the spatially 
extended coherence among the two-quasiparticle states. 
We call such spatially extended distribution of the two- 
quasiparticle states "the broad localization". 

In Figs E| and ^| the spatial structure of the two- 
quasiparticle states involving the 57/2 resonant state 
is also shown. Because the upper and lower compo- 
nents of the g r / 2 quasiparticle state are strongly confined 

(2) 

within the centrifugal barrier, the Fuv (r) concentrate 
only around the surface region. In general, we can expect 
the large collectivity by the broad localization with low-^ 
neutrons, and the effect is weaker in the two-quasiparticle 
states involving larger-^ states. 

The broad localization is realized owing to the broad- 
ening effect in the upper component of the quasiparticle 

(2) 

wave function. For comparison, the Ru V and the regions 
(2) (2) 

\Ruv — r\ < a U v without the broadening effect, that are 
calculated with the lower components obtained by HFB 
and the upper components obtained by the resonant BCS 
approximation, are shown in Fig ll5l If the broadening ef- 
fect is absent, the distributions concentrate only around 

(2\ 

the surface region with small a^v . In this situation, the 
collective motions among them can appear, however the 
low-frequency vibrational excitations have a similar be- 
havior with those in stable nuclei, and novel features of 
neutron drip line nuclei are suppressed. 



IV. HFB PLUS QRPA CALCULATION 

A. Formulation and inputs 

We consider the first 2 + states in neutron rich Ni iso- 
topes to examine the unique role of self-consistent pair- 
ing correlation in the neutron drip line nuclei. In Fig llSI 
the neutron single-particle energies obtained by HF cal- 
culation with Skyrme SLy4 force [25| are shown. The 
neutron drip line nucleus is 86 Ni within the HF calcula- 
tion. By taking into account pairing correlations, more 
neutrons can bound. However the predicted position of 
the neutron drip line depends on the treatment of pair- 
ing correlations, for example, the drip line nucleus is 88 Ni 
in Ref.|33j and 92 Ni in Rcf. 40J. The predicted drip line 
also depends on the effective interactions and the frame- 
works; namely relativistic or non-relativistic approaches 
(for example, see Ref . [40t l4l| ) . In the present study we 
consider up to 88 Ni, because our purpose is to investi- 
gate the qualitative aspects of vibrational excitations in 
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Mass Number 

FIG. 16: The neutron single-particle energies in neutron rich 
Ni isotopes obtained by HF with Skyrme SLy4 force. The 
approximate resonance energies with Rbox = 30 fm are also 
shown. The neutron Fermi energies in HFB are indicated by 
thick line. 



neutron drip line nuclei. 

We perform HFB plus QRPA calculation with Skyrme 
force. The HFB mean fields are determined self- 
consistently from an effective force and the residual in- 
teraction of the QRPA problem is derived from the same 
force. The QRPA problem is solved by the response func- 
tion method in coordinate space. A detailed account of 
the method can be found in Ref . |42l E3. l44 | . The residual 
interaction has an explicit momentum dependence, and 
these momentum dependence are explicitly treated in our 
QRPA calculations. Because we calculate only natural 
parity (non spin- flip) excitations, we drop the spin-spin 
part of the residual interaction. The residual Coulomb 
and residual spin-orbit interactions are also dropped. 

The ground states are obtained by Skyrme-HFB calcu- 
lation. The HFB equation is diagonalized on a Skyrme- 
HF basis represented in coordinate space with the box 
boundary condition of Rbox = 30 fm. Spherical sym- 
metry is imposed on the quasiparticle wave functions. 
The quasiparticle states with the energy below 50 MeV 
and the orbital angular momentum up to 7h are taken 
into account. The Skyrme SLy4 force is used for the 
HF mean-field. The density-dependent zero-range pair- 
ing interaction of Eq. Q with the density p(r) of protons 
plus neutrons is adopted for the pairing field. The pa- 



rameters are taken to be V m 



-555 MeV fm- J and 



that give the average neutron 
86 in 86 Ni. 



Pc = D.16(w_p(0)) fm- 3 
pairing gap A n w 12/ 



B. First 2 + states in Ni isotopes 

In FigH3 the B(E2, Of -> 2+) values and the exci- 
tation energies of the first 2 + states in neutron rich Ni 
isotopes obtained by HFB plus QRPA calculation are 
shown. These 2 + states are discrete solutions below the 
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FIG. 17: The B(E2, 0+ -> 2f) values and the excitation ener- 
gies of the first 2 + states in neutron rich Ni isotopes obtained 
by HFB plus QRPA, resonant BCS plus QRPA, and RPA 
with Skyrme SLy4 force. 
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FIG. 19: The isoscalar and charge quadrupole strength func- 
tions in 86 Ni obtained by HFB plus QRPA with Skyrme 
SLy4 force. The corresponding unperturbed strengths are 
also shown. A width parameter P = 0.05 MeV is introduced 
for the plots. The threshold energy E t h is indicated by arrow. 
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FIG. 18: The box size dependence of the B{E2) values ob- 
tained by HFB plus QRPA with Skyrme SLy4 force. 

threshold energies. The B(E2) values increase as ap- 
proaching the neutron drip line. The B(E2) values in 
the single-particle unit (=5 times the Weisskopf unit) 
are 7.8 in 86 Ni, and 11.7 in 88 Ni. To examine the con- 
vergence of the solutions, the box size dependence of the 
B(E2) values is shown in Fig^| Because the spatial 
distribution of the contributing two-quasiparticle states 
is localized by the pairing anti-halo effect, the coherence 
among them occurs within 20 fm as shown in Fig ll5l 
Therefore Rbox = 20 fm is enough for the convergence. 

In Fig lliJI the isoscalar and charge quadrupole response 
functions in 86 Ni obtained by HFB plus QRPA, and the 
corresponding unperturbed strength functions are shown. 



For the plots, a small width parameter of T = 0.05 MeV 
is introduced. All two-quasiparticle states are in the con- 
tinuum region, and the configurations below 6 MeV are 
of neutrons except for the proton (2p 3 / 2 )(l/7/2)~ 1 con- 
figuration at 4.9 MeV. The strong coherence among the 
neutrons, and the successive isoscalar coherent motion of 
protons and neutrons lead to the large B(E2) value. 

To examine the role of self-consistent pairing correla- 
tion, the results of the resonant BCS plus QRPA calcu- 
lation up to 86 Ni are also plotted in FigEI Although 
the resonant BCS approximation is not well-defined in 
neutron drip line region, the results are shown to empha- 
size the limited applicability of this approximation. The 
technical detail of our resonant BCS plus QRPA calcu- 
lation is explained in the previous report |37|. The box 
size of Rbox = 20 fm is used only for the resonant BCS 
approximation. The results have the smooth dependence 
on Rbox up to around 20 fm, and they are consistent 
with our preliminary results obtained by the resonant 
BCS plus QRPA calculation with the constant pairing 
gap approximation 37] ■ However we couldn't obtain the 
regular solution in 84 - 86 Ni with the larger Rb ox - In the 
present study, as an improvement, the density-dependent 
pairing interaction of Eq. is used. Because the cut-off 
energies of the particle model space are different in the 
resonant BCS calculation and the QRPA calculation, the 
different pairing strengths must be used for them. The 
pairing strength is changed for each nucleus to repro- 
duce the average neutron pairing gap obtained by HFB 
calculation. On the other hand, for the QRPA residual 
interaction, exactly the same interaction used in the HFB 
plus QRPA calculation is adopted. 
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FIG. 20: The isoscalar and charge quadrupole strength func- 
tions in 86 Ni obtained by RPA with Skyrme SLy4 force. The 
corresponding unperturbed strengths are also shown. The 
neutron lp-lh excitations; (a) 3si/ 2 — * resonant 0(3/2, (b) 
2^5/2 ~ * resonant d 3 / 2 , and (c) 2d 5 / 2 — > resonant gr/2, are in- 
dicated by vertical lines. The configuration, (x) 3s 1 / 2 — * (dis- 
cretized) non-resonant continuum d 5 /2, is also pointed out. A 
width parameter V — 0.05 MeV is introduced for the plots. 
The threshold energy E t h is indicated by arrow. 



Prom 80 Ni to 82 Ni, the results of the HFB plus QRPA 
and the resonant BCS plus QRPA have a similar behav- 
ior; namely the energies decrease and the B(E2) values 
increase. On the other hand, as approaching the neutron 
drip line, the B(E2) values obtained by the resonant BCS 
plus QRPA calculation are much smaller than those in 
the HFB plus QRPA calculation. The main qualitative 
difference between these calculations is the realization of 
the broad localization of the two-quasiparticle states as- 



2 3/2 



and d 



5/2 



sociated with the neutron low- lying Sj/2, 
resonant states in HFB. 

The results for 84 ' 86 Ni obtained by RPA calculation are 
also shown in Fig ll7l The B(E2) values in RPA is much 
smaller than those in the resonant BCS plus QRPA cal- 
culation, because the contributions of the hole-hole chan- 
nel can increase the collectivity of these excitations. As 

the lowest discrete so- 
The 



we have discussed in Subsec lHIBI 
lutions in RPA are single-particle like excitations, 
excitation energies are almost the same with the specific 
lp-lh energies (2d 5 / 2 — > 3sx/2 in 84 Ni, and 3siy 2 - ► res- 



onant d 3 / 2 in 86 Ni). In FigJ^Ul the isoscalar and charge 
quadrupole strengths in 86 Ni obtained by RPA calcula- 
tion, and the corresponding unperturbed strengths are 
shown. The neutron lp-lh excitations; (a) 3si/ 2 —> res- 



onant d. 



3/2, 



(b) 2d 



5/2 



resonant ^3/2, and (c) 2d, 



5/2 



— > resonant g-jn, are indicated. The excitation, (x) 
3s i/ 2 — + (discretized) non-resonant continuum d 5 /2, is 
also pointed out. 

The excitations to the resonant states, (a), (b) and (c), 
can generate some collectivity among the localized com- 
ponents (see Figsll2land ll3[) . However the B(E2) values 
of the first 2 + states in RPA change from 234 e 2 fm 4 to 
173 e 2 fm 4 in 84 Ni, from 83.0 e 2 fm 4 to 74.5 e 2 fm 4 in 86 Ni 
as increasing the box size from 20 fm to 30 fm. The HF 
calculation with small Rbox leads to the artificial local- 
ization of the spatially extended wave functions. Accord- 
ingly the RPA calculation overestimates the correlations 
among these lp-lh states. To obtain fully convergent 
solutions in RPA, much larger Rbox or continuum RPA 
calculation is required. The excitations to non-resonant 
continuum states like (x) can't produce the collectivity. 
The small peak related to (x) in the charge quadrupole 
strength is an artifact of the box boundary condition, 
and the strength decreases with larger Rb ox - 



V. CONCLUSION 

We have investigated the role of low-£ neutrons and 
pairing correlations for low-frequency vibrational exci- 
tations in neutron drip line nuclei. We have clarified 
the change of the spatial structure of quasiparticle wave 
functions induced by self-consistent pairing correlation; 
the pairing anti-halo effect in the lower component and 
the broadening effect in the upper component. We have 
found that the broad localization of two-quasiparticle 
states among low-^ neutrons is realized in consequence of 
the structure change in the quasiparticle wave functions. 
The broad localization can cause novel features of low- 
frequency vibrational excitations in nuclei close to the 
neutron drip line; as demonstrated by HFB plus QRPA 
calculation for the first 2 + states in neutron rich Ni iso- 
topes, it brings about the spatially extended coherence 
and the large transition strength. 
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